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Abstract 

In this paper we present the distinguished (d-) Riemannian geometry 
(in the sense of nonlinear connection, Cartan canonical linear connection, 
together with its d-torsions and d-curvatures) for a possible Lagrangian 
<— j , inspired by optics in non-uniform media. The corresponding equations 

O j of motion are also exposed, and some particular solutions are given. For 

instance, we obtain as geodesic trajectories some circular helices (depend- 
ing on an angular velocity cj), certain circles situated in some planes (ones 
are parallel with xOy, and other ones are orthogonal on xOy), or some 
straight lines which are parallel with the axis Oz. All these geometrical 
geodesies are very specific because they are completely determined by the 
non-constant index of refraction n(x). 

Mathematics Subject Classification (2010): 53C60, 53C80, 83C10. 
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1 Introduction 

In the geometrical optics [3], a central role is played by the Synge-Beil metric 

ID El El 

g a p{x,y) = <p aP (x) + i 2 y a yii, (1.1) 

where "f(x) > is a positive smooth function on the space-time M 4 , and <p a p(x) 
is a pseudo- Riemannian metric on M . One assumes that the four-dimensional 
manifold M A (which is connected and simply connected) is endowed with the 
local coordinates 

[x ) Q _Q 3 \x t, x , x , x ) ; 

for simplicity we use the system of units where the light velocity is c = 1. 
Obviously, the following rule holds good: y a — fa^y^- Because the components 
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of <p a s(x) are dimensionless, the same are fy a '- 

[fapi^)] = 1, hVa] = 1- 

In such a context, let us restrict our geometric-physical study to the Eu- 
clidean manifold £ 3 = (S 3 , Sij ) which has the local coordinates (a;) := (a: 1 )^^ • 

It follows that the corresponding tangent bundle TS 3 has the dimension equal 
to six, and its local coordinates are[§ 



(x,y) := {x\y l ) l=T - 



y spatial coordinates spatial direction J 



Let us introduce a metric on T£ 3 inspired by optics in a non-uniform medium 
(see formula (|1.1|> ): 

8ij(x,y) = <% +i 2 {x)y l y J , 

where 6 = (Sij) = diag (1, 1, 1) is the Euclidean metric, and yi = Si r y r . Usually, 
we have 7 2 (x) = n 2 (x) — 1, where n = n(x) is the refractive index of the non- 
uniform medium (see [H El [7]) ■ Using this spatial metric, below we will examine 
the special case of a possible anisotropic non-relativistic dynamical model, which 
is governed by the Lagrangian (in this model one considers that the particle has 
the mass m = 1) 



L(x,y) = -Qij(x,y)y l y J 



1 

^[% + 7 2 2/ l %]yV= (1.2) 



- |%yV + -VIM' 4 

where ||y|| 2 = (y 1 ) 2 + (y 2 ) 2 + (y 3 ) 2 = %yV . 



1 ■ ■ -y 2 

~ 2 w + \ 



Remark 1.1 Because the Euclidean metric 8^ is invariant with respect to the 
linear transformations of coordinates induced by the Lie group of orthogonal 
transformations 

0(3) = {AeM 3 (R)\A T -A = I 3 }, 

it immediately follows that the Lagrangian il.Sfy has a global geometrical char- 
acter with respect to these orthogonal transformations. 

Following as a pattern the geometrical ideas from Lagrangian geometry of 
tangent bundles [5] or jet bundles [5] , in what follows we construct the Riemann- 
Lagrange geometrical objects (the canonical nonlinear connection, the Cartan 
canonical linear connection, together with its d-torsions and d-curvatures) pro- 
duced by the Lagrangian (|1.2[) . 



§In this paper the Latin letters i, j, k, ... run from 1 to 3. The Einstein convention of 
summation is adopted all over this work. 
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2 Geometrical objects in the non-relativistic ex- 
tended dynamics 

The fundamental metrical distinguished tensor induced by the Lagrangian (|1.2[) 
is given by 



1 d 2 L 



a{x 1 y)8 ij + 2>y 2 {x)y l y 3 , 



2 difdyJ 

where <r{x,y) = (1/2) + J 2 {x)\\y\\ 2 > 0. 

Remark 2.1 Because the quadratic form q(£) — (yiyj)CS. 3 * s degenerate and 
has the signature (2,1,0), while the Euclidean metric 5(£) — Sij^ 1 ^ 3 is non- 
degenerate and has the signature (0,3,0), we easily deduce that the quadratic 
form <?(£) = gij{x 1 y)S ! l S > :l has the constant signature (0,3,0). It follows that it is 
invariant under orthogonal linear transformation of coordinates. Consequently, 
all the subsequent geometrical objects constructed in this paper will have the 
same form in any chart of coordinates induced by a linear transformation of 
coordinates produced by the orthogonal group 0(3). 

The inverse matrix g^ 1 = {g^) has the entries 

1 M 2 1 2 {x) 



f k (x,y) = 



;y J y k , 



a(x,y) a(x,y)-r(x,y)' 
where S jk = S jk andrfoj/) = (1/2) + 3 7 2 (a;)| |y| | 2 = a(x, y) + 2 7 2 (a;)||y|| 2 . 
Proposition 2.2 For the anisotropic Lagrangian il.ty) . the action 

E(x(t)) = [ L(x(t),y(t))dt, 

J a 

where y = dx/dt, produces on the tangent bundle TE 3 the canonical nonlinear 



connection N 



N) 



[Nj) , whose components are 
2 7 

—y % yj(i s y s )+ 

a 



7l|y| 



2 7 2 



5)(.i s y s ) + y l ij - i l yj y^jhsy") 



6 7 2 . 

— y l yjh s y s ) 

T 



(2.1) 



7 3 ibli 4 

2a 



1 3 3 

-Yyj - -y l ij - -$)h s y s )+ 

ST T J T~ ■* 



6t 2 

--\ T + 3 <7 )y l yj h s y s ) 



where j s = 9 7 jdx s and 7 l = <T r 7r 



7, 
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Proof. For the energy action functional E, the associated Euler-Lagrange equa- 
tions can be written in the equivalent form (see [3[2]) 

^f + 2G i (x k (t),y k (t))=0, V*=173, (2.2) 
where the local components 

a . def g ir r d 2 L . dL 
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dy r dx s dx r 



= -IMIV(7 s y s ) - -rIMIV - ^l|y|IV(7,y s ) 

a 4cr Zctt 

represent, from a geometrical point of view, a semispray on the tangent vector 
bundle TE 3 . The canonical nonlinear connection associated to this semispray 
has the components (see [5]) 

i dg dcr_ 

In conclusion, by direct computations, we find the expression (|2.1[) . ■ 
Remark 2.3 In an uniform medium with the constant refractive index 

n(x) =ti£ [1, oo), 

we have ~/ s — 0. Consequently, in this case we obtain G l — and Nj = 0. 



The nonlinear connection (|2.1|) produces the dual adapted bases of d- vector 
fields 

'^ = ^-^ 5 ^} CAr(TE8) (2 ' 3) 
and d-covector fields 

{dx i ; Sy l = dy* + N l r dx r ) C A-*(TE 3 ). (2.4) 

The naturalness of the geometrical adapted bases (|2.3j) and (|2.4j) is coming from 
the fact that, via a general transformation of coordinates, their elements trans- 
form as tensors on S 3 . Therefore, the description of all subsequent geometrical 
objects on the tangent space TE 3 (e.g., the Cartan canonical linear connection, 
its torsion and curvature) will be done in local adapted components. 
For instance, using the notations 

Nn := NT5 rj , N i0 := N„y r , N oj := N rj y r , N 00 := N^y*^, 

by direct local computations, we obtain the following geometrical result: 

Proposition 2.4 The Cartan canonical N-linear connection produced by the 
anisotropic Lagrangian il.ty) has the adapted local components 

CT(N) = (L%, C%) , 
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wh 



ere 



L \k = ~b (4^o + s i N oo ~ s 3k s ir N r0 ) + 1 Ml 2 (s Jk y - 5) lk - 5* fc7j ) + 



+7 {(N jk + N kj ) y l + (Nt - S tr N rk ) y 3 + (iVj - 5 ir N rj ) y k ) + 

2'Y 3 7/ i 

+2 (Yyjyk - y l y 3 i k - v l Vki 3 )} + [7 (ifeMw + ifc-fyo - tfjkiVoo) + 

err 

+lly|| 2 (S]k7 r y r - Vjik - Vkij) + 2 {y ] yki r y r - vnkWvW 2 - Vki d \\v\^) + 

+7 {(N jk + N kj ) \\y\\ 2 + {N k0 - N ok ) y 3 + (N j0 - N oj ) y k }] , 

-Y 2 2-y 4 
C )k = ^ ( S )yk + Slyj + 5 jk y l ) - -j- {\\y\\ 2 5jk + 2y 3 y k ) y\ 

Proof. The adapted components of the Cartan canonical connection are given 
by the general formulas (see [5]) 

T i def g ir ( Sg jr Sg kr Sg jk 



2 V Sx k 5xi 5x r 



n i def g ir ( dg jr dg kr dg jk \ _ g ir dg jr 



J]k 2 \dy k dyi dy r J 2 dy k 

Using the derivative operators (|2.3I) . the direct calculations lead us to the re- 
quired results. ■ 

The Cartan canonical iV-linear connection produced by the anisotropic La- 
grangian (|1.2p is characterized by three effective local torsion d-tensors, namely 

t dg SNl_6Nl pi defdN[_ , , 
'' 8x k 5x3 ' dy k kj ' jk ' 

and three effective local curvature d-tensors: 

jf,! 5L l - k 5L\, 

OX OX j j j 

pi def j k ~j ~j „ r 

jkl ~ ~~dy T ~ j7|fc jr kh 
Jkl ~ ~dy r ~ "dp 7 jk rl ~ j ' 1 rk ' 



where 



At the end of this Section, we would like to point out that in our anisotropic 
Lagrangian geometrical theory for non-relativistic extended dynamics, the clas- 
sical Riemannian Levi-Civita connection (which in our case would be attached 
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even to the Euclidean metric Sij) is replaced by Cartan canonical connection 
(associated to the perturbed Lagrangian (|1.2|l ). It is well known that in Finsler- 
Lagrange geometrical framework (see [5]) there are a lot of important linear 
distinguished (d-) connections (Cartan, Berwald, Chern-Rund or Hashiguchi, 
for instance). However, the use of the Cartan canonical connection is preferred 
because it is the single linear d-connection which is a metrical connection, like 
the Levi-Civita connection in the classical Riemannian framework. An impor- 
tant difference between these two connections is that the Cartan connection is 
only partial torsion- free because the Poisson brackets of the distinguished vector 
fields S/Sx l are not generally equal to zero. 

3 The equations of motion in the anisotropic 
non-relativist ic dynamics 

The Euler-Lagrange equations (|2.2p , for 



{y 1 ,y 2 ,y 3 ) = (V 1 ,V 2 ,V 3 ):=V, V 1 



dt 



v 2 = \V\ 2 = (V 1 ) 2 + (V 2 ) 2 + (V 3 ) 2 , 
lead us to the following anisotropic equations of motion for a particle: 

^ + 4 V (1 + 37V) , vsyi) _ jv* (31) 
where i E {1, 2, 3}. 

Remark 3.1 In a uniform medium with the constant refractive index 

nix) = no G [l,oo), the above anisotropic equations of motion simplify as 

dV i 

— = & V = (V 1 , V 2 , V 3 ) = constant 
dt v ' 

dx l 

— = V l & x{t) = (vh + xl, V 2 t + x 2 ,V 3 t + x 3 ) , 

where (xq,Xq,Xq) = constant. It follows that, in this case, the particles are 
moving only on straight lines. 

3.1 Particular solutions with cylindrical symmetry 

Let us investigate now the case of a non-uniform medium with cylindrical 
symmetry This means that we have n 2 = 1 + f 2 (p) ^7 = /(/?), where 
/ : (0, 00) — > (0, 00) is an arbitrary non-constant smooth function, and 

p 2 = (x 1 ) 2 + (x 2 ) 2 > 0. 



6 



In such a special context, let us search for solutions of (|3.ip in cylindrical 
coordinates: 

x 1 ^) = p{t) cos (f>{t), x 2 (t) = p{t) sin 4>(t) , x 3 (t) = ((t), 

where 6 € [0, 2ir] and (el. By direct computations, we deduce that the 
equations (|3.1[i rewrite as 



+ 



6(p) — 2 sin 6 (p<f>j — P cos 6 — p sin 6 Cj> \ 4 

4//V (l + 3/V) / a / .u 

l + 2fV)(l + 6W V COS ^ ( ^ -P***\M)) 



(l + 2/ 2 v 2 )(l + 6/ 2 ?; 2 ) 



(3.2) 



. cosfli = 0, 



1 + 2f 2 v 2 

sin 6 (p) + 2 cos f P<^ — p sin ^(/>^ + p cos </> 
4//V(l+3/V) 



(l + 2/ 2 w 2 )(l + 6/ 2 t; 2 ) 



sine/) (p) + pcos6 
sin (A = 0, 



where 



1 + 2fV 

4//V(l + 3/V) /.x 
^ (l + 2/ 2 w 2 )(l + 6/ 2 w 2 ) V^J 

v 2 = (p) 2 + P 2 W 2 +(C) 2 - 



0, 



(3.3) 



(3.4) 



In what follows, we will look for some particular solutions of the system of 
differential equations (13.21 - 13.41) . 

Case 1: Let us consider that p = and u ^ 0. We will look for solutions of 
the form 



c(C) = {p = constant, 6 = 6((), ( e 1) ^ c(C) = ( p = 0, 
These conditions imply 



d6 



C = i 



2 2 
V = p 



d6 



1. 



In this case the anisotropic equations of motion reduce to 

2 



< 



-poos < 



-psm < 



p sm ( 



p cos < 



d 2 6\ 


ffv i 


dC 2 ) 


1 + 2/ 2 w 2 


d 2 6\ 


ffv 4 


dC 2 ) 


1 + 2/ 2 w 2 



cos 6 = 



sin 6 = 0. 



(3.5) 
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Multiplying the first equation of (|3.5p by (—sin <j>), and the second by (cos</>), 
by summing we get the equation: 



'0 = 



(3.6) 



where w, O € R are arbitrary constants (here u> has the physical meaning of 
angular velocity). Consequently, the function (13.61) is a solution for the equations 
of motion (|3.5I) if and only if 



ffv 



'„,4 



+ 2/ 2 v 2 







(pcos</>) (w 2 ) + y 
(psin0) (w 2 ) + - ^ J 2 ^ 2 sin0 - 



<^> pW 



//V 



1 + 2/V 



= 0, 



(3.7) 



where v 2 = p 2 uj 2 + 1. 
1. If 2/ 2 - 1< and 



V + Pf e 



(2/ 2 - 1) 
4/ 



-,0 U(0,2/), 



then the solutions of the equation (I3.7[) are 



w = ±- 
P 



1 2f 2 -l±VA 



2f(2f + pf>) 



-1, 



where A = (l - 2/ 2 ) 2 + 4/ (2/ + p/') . 
2. If 2/ 2 - 1 > and 2/ + pf E (0, 2/) , then the solutions of ([3J]) are 



1 / 2/ 2 -l + 7A 

"° = :± pV wwtpY) ~ l 

In conclusion, the solutions of equations of motion (|3.5[) are 
c(t) = (p = p = constant, </> = cjgt + </>o, C = t) =>■ 

x x (t) = p o cos0(i), x 2 (<) = p sin<^(t), x 3 (<) = i. 

It follows that in this case the particles move on some circular helices. Note 
that these circular helices are some very specific trajectories (non-evident in 
advance) because the values of the angular velocity w are completely deter- 
mined by the initial function /(p). Solutions corresponding to the sign plus and 
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minus, which are intimately related to the orientation of rotation (right-handed 
and left-handed respectively). 

Case 2: Let it be </> = and v ^ 0. We will look for solutions of the form 



C (C) = (p = p(Q, <f> = constant, (eR)^ c(() = p 



dp 



o, c = i 



These conditions imply 



In this case the anisotropic equations of motion reduce to 



cos< 



sin < 



fp 



fp 



4//V(l + 3/V) 
(1 + 2f 2 v 2 ) (1 + 6f 2 v 2 ) 

ffv 4 , n 

~ 1 i Ot2 2 C0S $ = 

1 + 2j z v z 

Iff'v 2 (1 + 3/V) 
(1 + 2f 2 v 2 ) (1 + 6f 2 v 2 ) 

ffv* 



cos < 



sine 



(3- 



1 + 2f 2 v 2 



sin < 







The last equation of (|3.8p implies dp/dt^ = or /' = 0. In both situations the 
solutions of the equations of motion are 

c(t) = (p = p = constant, <j) = 4> £ [0, 2w], ( = f e 1) =^ 

=> = p cos0 o , x 2 (t) = p sin</> , x 3 (t)=t. 

Consequently, in this case the trajectories are the generators of the right 
circular cylinders 

(^) 2 + {x 2 f = pi 

where p > is a solution of the equation ,f'{p) = 0. 

Case 3: Let us consider that £ = and v ^ 0. We will look for solutions of 
the form 

c((f>) = (p = p((f)), <j> := t £ [0, 2tt], ( = constant) =>■ 
dp 



t{4>) = P 



l, c = o . 



These conditions imply 



dp 

*' + ' r 
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In this case the anisotropic equations of motion reduce to 
i 2 p\ „ . fdp\ ff'v i 

4//V(l + 3/V) ( fdpV . fdp\\ 

cost — -pant — =0 



(l + 2/ 2 w 2 )(l + 6/ 2 w 2 ) V \dt) r \dt 

. /d 2 p\ „ (dp\ . ff'v i . 
sin t —-5- + 2 cos 1 1 — I — p sin t 7o~o sm ^+ 



(3.9) 



di 2 J \dtj r 1 + 2/ 2 ?; 2 

4//v(i+3/v) f . ,/d P y w 



(l + 2/ 2 u 2 )(l + 6/ 2 w 2 ) I \dt J r \dt 

Multiplying the first equation of (|3.9[) by (— sin</>), and the second by (cos^>), 
by summing we get the equation: 

dp\ ( 2p//V(l + 3/V) \ = 
dt ) I (1 + 2/ 2 w 2 ) (1 + 6f 2 v 2 )J 

1. If dp/dt = 0, then the equations of motion (|3.9[) become 

p cos t H — — — cos f — 

P l + 2/ 2 w 2 ffV 

//V P 1 + 2/V 2 

^ sint+ 1 + 2/^2 smt = Q - 

^ 1 + 2/ V + //'p 3 = 0, (3.10) 
where v 2 — p 2 . The corresponding solutions are 

c(t) = (P = Po > 0, ^ = t G [0,2tt], C = Co G R) =► 

=>■ x 1 (<) = p cost, x 2 (i) = p sint, x 3 (t) = Co- 
It follows that in this case i/ie trajectories are circles situated in planes 
parallel with the plane xOy, having the centers situated on the 
axis Oz, and the radii p Q > 0, where p is a solution of (3.10\) . 

2. If dp/dt 7^ 0, then we get 

Zpff'v 2 (1 + 3/V) , 

ij i . l = n (311^ 

(l + 2/ 2 i; 2 )(l + 6/ 2 t> 2 ) V ; 

The preceding equation has solution if and only if 

2f+pf H"W(ITW)' 
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In this case, the solutions of the equation p. lip arc 



dp -4p- P ff-VN 2 

~di = \ wwtpT) — p :=F{ph (3 - 12) 

where A' = 4/ 4 + p 2 f 2 {f'f + 2pf 3 f. Imposing as a solution of (|3~i~2j) to 
be also a solution for the system (|3.9|) . we get 

FF ,_ p 2F 2 ff'(F 2 + p 2 ) 2 



p l + 2f 2 {F 2 +p 2 ) 

=>■ p = constant (contradiction!). 
So, the case dp/dt ^ does not provide us any solution. 

3.2 Particular solutions with spherical symmetry 

We investigate now the case of a non-uniform medium with spherical symmetry. 
This means that we have n 2 = 1 + / 2 (r) 7 = /(r), where / : (0, 00) — > (0, 00) 
is an arbitrary non-constant smooth function, and 

r 2 = |x| 2 = (x 1 ) 2 + (x 2 f + (x 3 ) 2 > 0. 

Example 3.2 In the study of mirages from heat above warm surfaces, one 
works with indices of refraction which have spherical symmetry (see JQWjj): 



n(x) = 1 + ee- r /a & f(r) = ^2ee-^l^ + e 2 e~ 2r2 /°\ 
where e > and a > are some given constants. 



In such a special context, let us search for solutions of (|3.ip in spherical 
coordinates: 

x x (t) = r(f) sin 0(f) cos 0(f), x 2 (t) = r(f) sin 0(f) sin 0(f), x 3 (f) = r(f) cos 0(f), 

where € [0,7r] and € [0, 27r]. By direct computations, we deduce that the 
equations (13.11) rewrite as 

sin cos (r) + r cos cos ( 0^ — r sin sin ^0^ + 

+2 cos 9 cos — 2 sin sin (rc/?J — 2r cos sin 



-r sm cos e 



4//V(l + 3/V 



(l + 2/ 2 w 2 )(l + 6/ 2 w 2 ) 
^sin cos (r) 2 + r cos cos — r sm $ sm ) " 



(3.13) 



1 + 2/ 2 u 2 



sin cos = 0, 
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sin 9 sin <f> (f) + r cos 9 sin <j> (() \ + r sin 9 cos 4> Cf) + 

j + 2r cos cos ^ft 
4//V (1 + 3/ V) 



+2 cos sin 



-r sin ft sin < 



2 sin cos ( f<; 

2 / . \ 2 X 



(1 + 2f 2 v 2 ) (1 + 6f 2 v 2 ) 



(3.14) 



^sin 9 sin (f) + r cos sin (j) y9j + r sin 9 cos 
//V 

"TTW sin0sin</, = 0, 

cos 9 (r) — r sin ^ft^ — 2 sin — r cos $ 
4//V(l + 3/V) 



^(^)) 



+ 



(1 + 2/V) (1 + 6f 2 v 2 ) 
ffv 4 



1 + 2f 2 v 2 



cos 6* (f) — r sin ( rftj ) — 
cos ft = 0, 



(3.15) 



where 



v 2 = {rf + r 2 (e) + r 2 sin 2 8 



Let us look for some particular solutions of the DEs system p. 131 - [3". 15|) . 
Case 1: Let it be f = and n ^ 0. We will look for solutions of the form 



c(<p) = (r ~ constant, 9 

^c((f>) = | r = 0. 
These conditions imply 



, <t>:=te [0,2tt]) 
d8 



1 



2 2 

f = r 



sm 



In this case the anisotropic equations of motion reduce to 
'd 2 9\ _ „ . fd8^ 



r cos ft cos t 



dt 2 



2r cos 9 sin t 



- C Sll) ft CDS / j | — 



1 - 



dt 
ffv 4 



1 + 2f 2 v 2 



sin ft cos t = 



,'d 2 ft\ /eZ0 
r cos sin i — — + 2r cos ft cos t — 



dt 2 



-rsin0sin£ 
'd 2 9 



r sin 9 



dt 2 



r cos 9 



1 



(ft 



(3.16) 



1 + 2f 2 v 2 
d9\ 2 , //V 



1 + 2/V 



sin ft sin f = 



cos0 = 0. 
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Multiplying the first equation of (|3.16p by (— sint), and the second by (cost), 
by summing we get the equation: 



2r cos 8 



dd\ 
dt) 



= 0. 



1. If cosd = 0-^8 = 7r/2, then the equations of motion (|3 . 16[) become 

ffv 4 



r sin 8 cos t + 
r sin 8 sint + 

1 



1 + 2f 2 v 2 

ffv 4 
1 + 2f 2 v 2 

ffr 3 



sin 8 cos t = 



sin 8 sint = 0, 



1 + 2/ 2 r 2 



(3.17) 



where v 2 



It follows that in this case the solutions of the anisotropic 



equations of motion are 



c(t) 



ro = constant, 8 = — , 



.x 1 (t) = r cost, x 2 (t) — r sint, a; 3 (t) = 0, 

where ro > is a solution of the equation (|3.17[) . It follows that in this case 
the particles move on the circles situated in the plane xOy, which 
have the centers in origin and the radii equal to the roots of 
the equation 1T\) . We recall that in the case of cylinder solutions 
we obtained some circles situated in planes parallel with the plane xOy, 
having the centers situated on the axis Oz, and the radii p > 0, where 
p is a solution of the same equation. 

2. If we have d6/dt — 0, then the anisotropic equations of motion become 

ff'v 4 



r sin 8 cos t + 
r sin 8 sint + 



ffv 4 



2f 2 v 2 



1 + 2f 2 v 2 

ffv 4 
+ 1 + 2f 2 v 2 

cos 8 = 0. 



sin 8 cos t = 
sin 8 sin t = 



This system implies 8 = tt/2 (we recover then the above Case 1.) or 
/' = 0. So, the new corresponding solutions are 

c(t) - (r - r > 0, 8 = 8 a G [0, tt]\{^/2}, 4> = t)=> 

=>■ a: 1 (t) = (r sin^o) cost, x 2 (t) — (r sin O ) sint, x 3 (t) = r cos#o- 
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These trajectories are circles situated in planes parallel with the 
plane xOy, having the centers situated on the axis Oz, and the 
radii r s'm9 , where r > is a solution of the equation f'(r) = 0. 
By comparing, we recall that in the case of cylinder solutions we obtained 
the generators of the right circular cylinders 

(X 1 ) 2 + (* 2 ) 2 = pi 

where p > is a solution of the same equation. 

Case 2: Let us consider that 9 = and ti^O. We will look for solutions of 
the form 

c{<j>) = (r = r(4>), 9 = constant, (j> := t e [0, 2tt]) 

=*m= (r= |^ = ' '<t>= 1 
These conditions imply 

" 2 = © 2 + r28i " 2 *' 

In this case the anisotropic equations of motion reduce to 

f d 2 r\ ( dr\ 

sin 9 cos t -— — 2 sin 9 sin t \ — — r sin 9 cos t+ 
\dt 2 J \dt J 

4//V(l + 3/V) ( . a (drV . a ..(dr 

+-r. ,-. — nx r- „ — nx sm " cos M ~r \ — r sin 9 sin t — 

(l + 2f 2 v 2 )(l + 6f 2 v 2 ) y \dt) \dt 

ff'v 4 

■ sin 9 cost = 0, 



1 + 2/V 

sin 6 sin t [ \ + 2 sin 9 cos t I — ) — r sin 6 sin t+ 

4//V(l + 3/V) / . /dr\ 2 /dr 
1 sin 9 sin £ — + r sin cos i 1 



(i + 2pv 2 ) (i + epv 2 ) y \dt J ' v dt 

ff'v 4 



1 + 2f 2 v 2 



sin 9 sin t = 0, 



Q (d 2 r\ , 4//V(l + 3/V) /dr 
C ° S * (d? J + (1 + 2/^) (1 + 6,/^) C ° S 9 U 

f fV 

" I I 0/2 2 COSg = 0. 

1 + 2/ 2 v J 



2 
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1. If cos = <^=> = it/2, then the equations of motion become 
'd 2 r\ „ . . fdr\ . ff'v 4 



cost 



dt 2 J 



- 2sini 



r cos t 



1 + 2/V 



cos t+ 



+ 



Affv 2 (1 + 3/V) 
(l + 2/ 2 v 2 )(l + 6/ 2 v 2 )' 

2 



(3.18) 



cosi 



— r sin i 



dr\ 
) 



sin £ I — -7T ) + 2 cos t I ^ 



dt 2 



di 



r sin i 



11=0, 

ff'v* 
l + 2f 2 v 2 



sint+ 



+ 



4//V (1 + 3/V) 
(1 + 2f 2 v 2 ) {1 + 6f 2 v 2 )' 

2 



(3.19) 



dr- 



ill 



sin i — + r cos i — = 



dr 



where 



df 



dr 

* 1 1 



Multiplying the first equation by (— sini), and the second by (cost), by 
summing we get the equation: 



2r//V(l + 3/V) 
(l + 2,fV)(l + 6/ 2 v 2 ) 



0. 



(3.20) 



(a) If dr/dt = 0, then the equations of motion become 



ffv 4 

r cos t + - „„ „ cos t = 



r sin i 



4^ r + 



1 + 2/ 2 w 2 

ff'v* 
1 + 2/ 2 w 2 

ffv 4 



<S=> 



sin t = 



<^> 



1 + 2/ 2 w 2 
<=> 1 + 2/V + ffr 3 = 0, 
where v 2 = r 2 . The corresponding solutions are 



(3.21) 



c(t) = (r = r > 0, 6 = |, tf> = i) 



a; 1 (t) = r cos t, x 2 (t) = r sini, x 3 (t) = 0. 

It follows that in this case the particles move on circles situated in 
the plane xOy, having the centers in origin, and the radii 
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p > 0, where p is a solution of the equation ^3.21\) . Note 
that these trajectories also appear in the case of cylindric solutions 
(see again the preceding Case 3. // Subcase 1.). 

(b) If dr/dt 0, then we get 

2rff'v 2 (1 + 3/V) 

1 H — i '— = 

(1 + 2/V) (1 + 6/V) 

The preceding equation has solution if and only if 

/ 1 + 4r 2 f 2 

and the solutions are 



dr ^ -Ap-rff'-VN 2 

dt =± i 6/3(2/ + r/0 r ^ 

where A' = 4/ 4 + r 2 / 2 (/') + 2rf 3 f. Imposing as a solution of 
(|3.22[) to be also a solution for equations fj3. 18[) and (|3.19[) . we get 

p w 2^ 2 //'(^ 2 + r 2 ) 2 

i* i* — r — 



r l + 2/ 2 (F 2 +r 2 ) 

=>■ r = constant (contradiction!). 
So, the subcase dr/dt 7^ does not provide us any solution. 

2. If cos6> 7^ and sin^ 7^ 0, then the equations of motion reduce to the 
equations (|37T8j) . (|3~19f and 

, 4//V (1 + 3/V) fdr\ 2 



dt 2 J (l + 2/V)(l + 6/V) \dtj ( 

ff'v 4 
— = 0. 

1 + 2/V 

Using the equation (I3.23|) . we deduce that the equations f|3 . 1 8[) and (|3.19l) 
simplify as 

..fdr\ . 4r//V (1 + 3/V) , Jdr\ 

— 2sinr — — rcost — „„ „; ; ;^tV smt \ — =0, 

\dt J (l + 2/V)(l + 6/V) \dtj 

Jdr\ . ,, 4r//V (1 + 3/V) ( dr\ 

2 cos t\ — — r sin t + „ „ , . tttV cos t — = 0. 

\dt J (1 + 2/V) (1 + 6/V) V d V 

Multiplying the first equation from above by (cost), and the second by 
(sin t) , we deduce that this system of equations has no solution. 
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3. If sm# = 0, then the equations of motion are equivalent with differential 
equation f|3 . 23[) . where 



v 2 



To integrate the equation (|3.23[) seems to be very difficult but it is obvious 
that the solutions of the equations of motion are segments of the axis 
Oz, whose lengths are determined by the images of the solutions r(t) of 
the equation (I3.23[) . 

Case 3: Let us consider that <f> — and v ^ 0. We will look for solutions 
of the form 

c(9) = (r = r(9), 6 := t G [0,7r], <f> = constant) 

^c(6)= (r=^, 8 = 1, 4> = 
These conditions imply 

In this case the anisotropic equations of motion reduce to 

, fd 2 r\ „ , fdr\ 

sin t cos <p I 1+2 cos t cos <p I — I — r sin t cos < 



4//V(l + 3/V) ( ,(dr\\ (dr\ 

1 sm t cos © — + r cos t cos © — 



2 

(1 + 2f 2 v 2 ) (1 + 6/ 2 w 2 ) 1*"" "^^\dt) 1 ' " ^ v ^ y 



1 + 2/ 2 i; 2 



sin £ cos 6 = 0, 



• , ^ 2 r\ „ • , (dr\ . . , 
sin t sm <p I -^j 1+2 cos £ sm <p I — I — r sm t sm 0+ 

4//V(l + 3/V) f ,(dr\\ t . ,(dr\ 

1 sm t sm © — + r cos c sm © I — 



(l + 2/ 2 w 2 )(l + 6/ 2 w 2 ) \ T \dtJ \dt J 

ff'v A 



1 + 2/ 2 v 2 



r 



sin t sin = 0, 



dr \ 



cost |_l_ 2s intl-]-rc OS t- 



4/jV(l + 3/V) / /dr\ 2 . /dr 
1 cos i — — r sm t 1 



(l + 2/ 2 « 2 )(l + 6/ 2 « 2 ) ^ UB6 V^J """"1* 

//v 



1 + 2/ 2 w 2 



cost = 0. 
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Multiplying the first equation by (cos </> cost), the second by (sin </> cost), and 
the third by (— sint), by summing we get again the equation (|3.20[) . But, this 
case was previously treated, and the corresponding solutions are 



=> a; 1 (t) = (ro cos</> ) sinf, x 2 (t) = (r s'm(j) ) smt, x 3 (t) = r cost, 

where r > is a solution of the equation p. 211) . It follows that in this case the 
trajectories are the following circles: 



4 Conclusion 

At the end of this paper, we consider it is important to underline some geo- 
metrical differences between the classical non-relativistic dynamics (determined 
only by the Euclidian metric Sij) and the present anisotropic non-relativistic 
extended dynamics (determined by the perturbed Lagrangian (jl.2p ). For in- 
stance, we recall that in the first case we have only a Riemannian curvature 
(determined by the Levi-Civita connection), while in the second case we have 
more Lagrangian curvatures (determined by the more complicated Cartan con- 
nection). Moreover, in the first case we live in a flat space (its Riemannian 
curvature is zero), while in the second case we work in a curved spaces (its 
Lagrangian curvatures being non-zero) . At the same time, note that the corre- 
sponding geodesies are in these different cases: 

1. any straight lines in space, in the first case; 

2. only certain circular helices or circles situated in some specific planes, 
together with certain specific straight lines which are parallel with the 
axis Oz, in the second case. It is also important to note that, in this 
second anisotropic situation, the preceding geometrical geodesies do not 
represent the general solutions of the geodesic equations, but only some 
particular cases of them. To find the general solution of these geodesic 
equations can be considered as an open problem for a research study in 
this geometric-physical domain. 
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